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ABSTRACT 

Recent X-ray observations have revealed a universal temperature profile of the in- 
tracluster gas of non-cooling flow clusters which is flat for r < 0. 27-180 ■ Numerical 
simulations, however, obtain a steeper temperature profile in the inner region. We 
study the effect of thermal conduction on the intracluster gas in non-cooling flow clus- 
ters in light of these observations, using the steep temperature profiles obtained by 
authors of numerical simulations. We find that given 10^*^ yr for the intracluster gas to 
evolve, thermal conduction should be suppressed from the classical value by a factor 
^ 10^'^ in order to explain the observations. 



Key words: Cosmology: Theory — Galaxies: Intergalactic Medium- 
ters : general — X-rays: Galaxies: Clusters 



Galaxies : clus- 



1 INTRODUCTION 

The role of thermal conduction in the intracluster medium 
has long remained uncertain. For clusters of galaxies which 
show signs of cooling flow, it has been invoked in the past 
(Bertschinger & Meiksin 1986; David, Hughes & Tucker 
1992) to allow multiphase cooling flows to exist, although 
the degree of conduction was assumed to be suppressed be- 
low the classical value (Binney & Cowie 1981), possibly as a 
result of tangled magnetic fields (e.g. Fabian 1994). In view 
of the recent observations of cooling flows, several authors 
have revived the idea of thermal conduction, albeit with a 
suppression factor of only a few (Narayan & Medvedev 2001; 
Voigt et al. 2002). 

Concurrently, there have been two important develop- 
ments in regard to the temperature profile of non-cooling 
flow clusters. On one hand, observations with BeppoSAX 
have revealed a universal temperature profile, normalised 
by the emission weighted temperature of the whole cluster, 
which shows flattening within r < 0.2 nso (Molendi & de 
Grandi 2001). On the other hand, numerical simulations ob- 
tain a universal temperature profile that is somewhat steeper 
in this inner region (Loken et al. 2002) . One obvious physi- 
cal mechanism that can flatten the temperature profile and 
that has not been included in the simulation is thermal con- 
duction (Suginohara & Ostriker 1998; Molendi & de Grandi 
2002; Loken et al. 2002). 

Recently Loeb (2002) have argued on the basis of com- 
parison of conduction timescale with the age of clusters that 
thermal conduction should be suppressed at least by a factor 
~ 0.15 in order not to allow drastic cooling of the gas in the 



cluster cores (of non-cooling flow clusters). Ettori & Fabian 
(2000) argued that the observed temperature jumps in the 
cluster A2142 require the thermal conduction across cold 
fronts should be suppressed by a factor larger than ~ 100. 
Vikhlinin et al. (2001) also argued in favour of a suppression 
factor of ~ 100. Even in the bulk of the gas, Markevitch et 
al. (2002) pointed out that a suppression factor larger than 
~ 10 is needed to explain the existence of small scale inho- 
mogeneities in the temperature profile of some clusters. 

The physical mechanism for and the degree of suppres- 
sion are still unclear though. The proposed mechanisms in- 
clude tangled magnetic field (Tribble 1989), plasma insta- 
bilities (Pistinner, Levinson & Eichler 1996, and references 
therein). In the case of suppression by fluctuations in the 
magnetic field, Chandran & Cowley (1998) found a sup- 



pression factor of order 10 



-10" 



Pistinner, Levinson & 



Eichler (1996) also advocated an inhibition factor of order 
~ 10"'^ as a result of electromagnetic instabilities. Recently 
Narayan & Medvedev (2001) extended the analysis of Chan- 
dran & Cowley (1998) adding more wavevectors to the fluc- 
tuation spectrum of magnetic fields and found a suppression 
factor of order ~ 5. 

In this paper, we calculate the time evolution of the 
temperature profile of the intracluster gas in non-cooling 
flow clusters in the presence of thermal conduction, with 
various degrees of suppression, assuming quasi-hydrostatic 
equilibrium, and constrain the suppression factor by com- 
paring the resulting temperature profile with the BeppoSAX 
data. 

We first discuss the temperature and density profiles of 
the intracluster medium and the role of thermal conduction 
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in §2. Wc then compare the conduction timescale with the 
age of clusters in §3. We present the main result of our work 
in §4. Throughout the paper, we use the ACDM cosmological 
model, with Qm = 0.3, = 0.7 and h = 0.65. 



2 HEAT CONDUCTION AND THE 
INTRACLUSTER GAS 

Heat is conducted along the gradient of the electron tem- 
perature in a plasma. Thermal conduction is termed unsat- 
urated if the mean free path of electrons is smaller than the 
scale length of the temperature gradient. It is said to 'satu- 
rate' otherwise. In the case of unsaturated thermal conduc- 
tion, the heat flux is given by, 



Q = -KVTe , 



(1) 



where the thermal conductivity of a hydrogen plasma is 
given by (Spitzer 1962), 



Ksp = 1.8 X 10"^Ti^^^(ln A)r^ergs"^ cm"^ K"^ , 



(2) 



where the Coulomb logarithm is (Sarazin 1986) (with rie as 
the electron density). 



(InA)e = 37.8 -h In 



( 108 k) ( 10-3 cm-s) 



-1/2 



(3) 



In the case of intracluster gas, we will find that the scale 
length of the temperature gradient for the relevant temper- 
ature profile is larger than the mean free path of electrons 
(~ 23kpc(r/10*/S')2(ne/10-^cm-3)-i (Sarazin 1986)) in 
general. This means that the above mentioned expression 
for unsaturated heat flux is relevant here. In fact, with the 
suppression of thermal conduction, as is discussed below, the 
mean free path is smaller than the above value, bolstering 
the case for unsaturated heating due to thermal conduction. 

Thermal conduction will tend to transport heat from 
the hotter (usually the inner) regions of the intracluster gas 
to the colder (usually outer) parts. If the speed of the en- 
suing gas flow is assumed to be slow, or, equivalently, if the 
gas density profile is assumed to be in an approximate sta- 
tionary state, then one can write, 

3pg{r)kB dT{r,t) 
2' 



dt 



= -V-Q = V-(KVT(r,t)), 



(4) 



where pg{r),kB, fJ^mp are the gas density, Boltzmann con- 
stant, mean molecular weight and the proton mass respec- 
tively. 

2.1 Temperature profile 

Recently, Molendi & de Grandi (2002) have used the ob- 
servations of ten non-cooling flow clusters with BeppoSAX 
to derive a universal temperature profile of the intraclus- 
ter gas. They have normalised the profiles by the emission 
weighted temperature of the cluster, and plotted the profiles 
as a function of r/r iso, where riso is the radius where the 
mean overdensity of the cluster is 180 times that of the ambi- 
ent density. The profiles are characterised by an isothermal 
core that extends to r ~ 0.2ri8(), beyond which the temper- 
ature decUnes in the outer region. 

Molendi & de Grandi (2002) have compared their ob- 
served temperature profile with those predicted by several 



numerical simulations. Most of the predicted profiles are fiat- 
ter than the observed profile in the outer region (Evrard et 
al. 1996; Eke et al. 1998; Bialek et al. 2001). The high reso- 
lution simulation by Frenk et al. (1999) did predict a tem- 
perature profile that matches well the observed profile in 
the outer regions. Recently, Loken et al. (2002) have per- 
formed a numerical simulation with the highest resolution 
to date, and determined a universal temperature profile that 
matches the observed profile of Molendi &z de Grandi (2002) 
very well in the outer region (r > 0.2ri8()). They found a fit 
to their universal profile that is given by (with (T) as the 
emission weighted temperature) , which we use as the initial 
temperature profile. 



Ti{r) = 



1.3S{Ti] 



(1. + ^ 



1.5r M .1 



(5) 



2.2 Gas and dark matter density profile 

We assume the gas to be in a quasi-hydrostatic equilibrium 

state as the temperature profile of the gas evolves. The initial 
density profile was calculated using the background dark 
matter density profile of Navarro, Frenk & White (1996) 
assuming hydrostatic equilibrium. The dark matter density 
profile is then given by. 



pdm{r) 



(6) 



where ps is a normalizing density parameter. The charac- 
teristic radius re is related to the virial radius ryir by the 
'concentration parameter' (c), as 

c=^. (7) 

The total mass of the cluster is assumed to be the mass 
inside its virial radius. The virial radius is calculated in the 
spherical collapse model to be, 



My. 



{47r/3)Ac{z)pc{z) 



1/3 



(8) 



where Ac(z) is the spherical overdensity of the virialized 
halo within rvir at z in the units of the critical density of the 
universe pc{z)- For our adopted cosmological model Sim = 
0.3 and Ha = 0.7, the value of Ac{z) is 100 (Komatsu & 
Seljak 2002). The total dark matter mass within a radius r 



M(< r) = A-Kpsr^sm{r/rs) , 



(9) 



where, m.{x) = ln(l -|- a;) — 

We first determine the initial gas density profile cor- 
responding to the initial temperature profile as discussed 
above, assuming hydrostatic equilibrium. 



d_ 
dr ' 



\ urriri ) 



GM{< r) 



(10) 



where pi is the initial gas mass density and M(< r) is the 
total mass within a radius r as given by equation (9). We 
normalize the gas density by requiring the total gas mass 
(within T\,ir) to be a fraction fg = Q.h/Q.m = 0.15 of the total 
mass of the cluster, where we have used Q,b = 0.02/1^^ as 
constrained by primordial nucleosynthesis (Buries & Tytler 
1998). 
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3 COMPARISON OF TIME SCALES 

One can define a conduction time scale for a splierically 
symmetric cluster as (Sarazin 1986; Suginohara & Ostriker 
1998; Loeb 2002), 

3 neksT 3 nefcsT 



2|V-(kVT)| 2|if(r-2Kf)| 



(11) 



This can be compared with the age of the cluster, which 
we assume to be 10^° yr (as in, e.g., Loeb 2002). Here the 
conduction coefficient k is supposed to be a fraction r? of its 
classical value Ksp (cquation(2)). In Figure 1 wo show (by 
solid lines) this time scale for the case of r; = 0.1 for clusters 
of masses 5 x 10^'*, 10^^, 2 x 10^^ M©, as a function of r/rvir, 
and where the dotted line denotes the age of the clusters. 
We have used the initial temperature and density profiles 
as discussed in the preceding section. The discontinuity in 
the solid lines occur as a result of the modulus in the previ- 
ous equation. It is seen from the comparison of time scales 
with 10^° yr that for 77 = 0.1, thermal conduction plays an 
important role in determining the temperature profile for 
r/vvir < 0.1-0.2 for the above range of cluster masses. It 
could be argued that for the cluster temperature profiles to 
flatten below r/r^ir < 0.2 in ~ 10^° yr, the suppression 
factor for thermal conduction should be < 0.1. 

Loeb (2002) argued on the basis of a similar 
comparison — constraining tcond/lO^^yr < 0.3 at rmo that 
7] < 0.15(T/10 keV)~^^^ in order for the cluster cores not to 
have cooled substantially over the last Hubble time. We note 
here that the emission weighted temperature of clusters is 
related to the total mass as M oc T^/^ (Evrard et al. 1996), 
so that the limit on 77 scales as M^^ . 

Although the comparison of timescales is a good indica- 
tor, a much better fimit on the suppression factor is expected 
from the calculation of the time evolution the temperature 
profile as a result of thermal conduction. For example, from 
Figure 1 it is seen that for rj = 0.1, say for a cluster with 
mass 5 X 10" Mq, tcond = 10^° yr at r/vyir ~ 0.1. This only 
means that the temperature would drop at least by a factor 
(1 — e~^) at that radius. This would, however, mean that the 
temperature profile would be flat within a radius not simply 
equal to r/r vir ~ 0.1 but considerably larger than this. In 
other words, in reality the suppression factor could be much 
smaller than that obtained by a simple comparison of time 
scales. 



4 EVOLUTION OF TEMPERATURE PROFILE 

Wo next determine the evolution of the temperature profile 
using equation (4) with varying degrees of suppression of the 
value of K. For numerical integration of the equation with 
Crank-Nicholson method, it is helpful to change the variable 
from T to T^/^ because k oc T^/^ (equation(2)). The most 
natural boundary condition is to keep the temperature at 
Tvir a constant, where the effect of thermal conduction would 
be minimum. We have therefore integrated equation(4), with 
the boundary condition that the temperature at Vvir is a 
constant. 

To facilitate comparison with the BeppoSAX tempera- 
ture profile, which is normalised by the emission weighted 
temperature, we have calculated the emission weighted tem- 
perature for the time-evolved profiles, within r'25oo as. 
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Figure 1. Conduction time scale is compared to the age of the 
cluster of lO-'" yr, for three clusters of masses 5 X 10^^, 10^ ', 2 X 
10^^ Mq for a suppression factor of 17 = 0.1. The regions of heat- 
ing and cooling of the gas are marked. The dotted line indicates 
the timescale of 10^" yr that relates to the cluster age. 



{T) = 



/;^™ 47rr2n(r)2eo.5-ior(r)dr 



(12) 



where n represents gas particle density and eo.s-io denotes 
the emissivity for a metallicity of Z/Zq = 0.3 relevant 
for the 0.5-10 keV band, which we calculated using the 
Raymond-Smith code. We note that the initial profile of 
Lokcn ct al. (2002), the emission weighted temperature was 
calculated within a radius of lh~^ Mpc. For this profile, 
we have adopted the empirical scaling relation of mass and 
emission weighted temperature of Finogucnov ct al. (2001). 

We plot in Figure 2 the results of numerically inte- 
grating equation (4) beginning with the initial temperature 
profile of Loken et al. (2002) and the corresponding density 
profile, computed on the basis of hydrostatic equilibrium. 
The left panels show the temperature profile at the end of 
10^° yr for 77 = 10"\ 10"^ 10"^ (dotted lines) and the right 
panels show the same for temperature normalized by the 
corresponding emission weighted temperature. We plot the 
profiles for M = 2 X 10^^ and 5 x lO^"* M©. The mass range 
is chosen to probe clusters with (T) in the range of ~ 4-10 
keV. 

The figures show that given an age of clusters of order 

~ 10^" yr, a suppression factor for the thermal conduction 
that is r; ^ 10^'^ produces temperature profiles that are flat 
for a large fraction of the inner region and that are incon- 
sistent with BeppoSAX observations. The discrepancies for 
a given suppression factor become larger for larger clusters, 
as expected since thermal conduction depends on T"'^^. It 
is, however, seen that a value of r; ~ 10~^ results in a tem- 
perature profile that is consistent with the BeppoSAX data, 
almost independent of the cluster mass. 
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Figure 2. Evolution of the temperature profile is shown 
for different suppression factors of thermal conduction r) = 
10~^, 10~^, 10~^ (dotted lines), beginning with the initial pro- 
file shown by the solid line, given 10^" yr of evolution, for 
M = 2 X 10^^ M0 (upper panels) and M = 5 X 10^* Mq (lower 
panels). The left panels plot temperature against r/nso and the 
right panels plot the temperature normalized by the correspond- 
ing emission weighted temperature. 



5 DISCUSSION 

5.1 Luminosity evolution 

We can also attempt to constrain the suppression factor 
from the fact that a substantial evolution in the temperature 
profile of a cluster over cosmological time scale would also 
change its X-ray luminosity by a large amount in that time 
scale. Recent observations, however, find that the X-ray lu- 
minosity of clusters scale as Lx oc (l-|-z)^'^ which is expected 
from clusters collapsing at various cosmological epochs (see, 
e.g., Arnaud et al. (2002)). As oc p'^R^T^^^, then with 
T oc {pR^)/R and p oc (1 -I- z)^, one finds a scaling of Lx 
at a given temperature with redshift as (1 + z)^'''. Viklilinin 
et al. (2002) found with CHANDRA that for a sample of 
22 clusters with z > 0.4, the X-ray luminosity for a fixed 
temperature scaled approximately as (1 -f z)^'^^"'^. If ther- 
mal conduction changed the temperature substantially then 
these observations could be at variance with these simple 
cosmological scalings. 

For the purpose of illustrating the point, consider a clus- 
ter of a given mass collapsing at z — 0.5, with an initial pro- 
file similar to the one we considered in the previous section 
(and a corresponding density profile from the consideration 
of hydrostatic equilibrium, taking into account the changes 
in the dark matter potential and in the fraction fli/flrn, 
appropriate for a collapse redshift of z = 0.5). If wc now 
let this cluster passively evolve with thermal conduction at 
work, then wc could track its evolution in the Lx —T space. 
Vikhlinin et al. (2002) have shown that the data correlate 
well if the emission weighted temperature is plotted against 
Lx X {I + z)^''' (their Figure 2b). In Figure 3, wc have plot- 
ted the evolutionary track of two representative clusters (of 
masses 10^^ and 2 x 10^® M©) from z = 0.5 to z = 0.0, for 
the same parameters. We have calculated L^ within a radius 



of 2 Mpc and the emission weighted temperature within a 
shell of 0.5-1 Mpc, as Vikhlinin et al. (2002) have done. We 
also show the data points from Vikhlinin et al. (2002) as well 
as their best fit by the dotted line. 

The exact calibration of Lx and emission weighted tem- 
perature is not important here, as the initial temperature 
(density) profile has been chosen for the purpose of illustra- 
tion, not taking into account how the 'universal' tempera^ 
ture profile could evolve in redshift. The important points to 
note in this figure arc the extent and direction of the evolu- 
tionary tracks. The tracks show how thermal evolution could 
introduce scatter in this plot. As Voit (2000) has argued the 
difference between the coUpaso redshift of a cluster and the 
redshift at which it is observed can introduce scatter in the 
plots relating the X-ray parameters. 

A small value of thermal conduction coefficient would 
only move it along the Lx{l + z)^'^ axis, mainly because of 
the redshift factor. Although thermal conduction only con- 
ducts heat from the inner to the outer part of the cluster, 
and on average the gas temperature docs not change, but 
since X-ray emission depends on the density, temperatures 
determined from X-ray emission would be lowered by ther- 
mal conduction. A large value of the cocfRciont would result 
in a track in which both Lx and temperature are changed 
with time. The data at present still have a large scatter and 
it is not possible to constrain the value of thermal conduc- 
tion from this, but it might be possible to do so in the near 
future, if the uncertainties in the temperature measurements 
decrease the errorbars in the T direction. 

Similar evolutionary tracks could also be used for other 
X-ray parameters, e.g., in the M — T space. The scatter in 
the data for total mass and the emission weighted tempera- 
ture, especially for high redshift clusters, such as in Schindler 
(1999), however is too large to allow any useful constraint on 
thermal conduction. As Loeb (2002) has pointed out a sub- 
stantial evolution in the temperature profile would also in- 
troduce a large scatter in the measurements of cluster abun- 
dance. While these constraints would become better with 
more data in the future, we have obtained a much stronger 
limit from the measurement of the temperature profile as 
described in the previous section. 



5.2 Limitations 

Finally, we wish to reiterate the basic assumptions made in 
our calculations to remind ourselves of the limitations of the 
results. The gas was assumed to be in an approximate steady 
state condition. We have also done our calculations relaxing 

this assumption somewhat, by requiring the gas density pro- 
file to be in hydrostatic equilibrium at every time step. This 
did not alter the results much, and we believe that our re- 
sults should be robust in this regard, as long as there is no 
large scale gas flow. It should also be pointed out that we 
have assumed a single value of the suppression factor for 
the whole of the intraclustcr gas, independent of density or 
temperature. 

We assumed the intracluster gas to evolve passively in 
time, without any cooling and heating sources acting in it. It 
is possible that heating sources, such as active galactic rmclci 
(Bohringer et al. 2002), which could be more prevalent in the 
inner region, could compensate the cooling due to thermal 
conduction in this region to some extent. When radiative 
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Figure 3. Evolution of a cluster collapsing at z = 0.5 in the space 
oi Lx X (l+z)~^-^ and emission weighted temperature is shown for 
two clusters with M = lO^^ and 2 x lO^^ Solid, dotted and 
dashed lines show the evolution for ri = 10~^, 10~^, 10~^. Filled 
circles along the tracks indicate the luminosity and temperature 
at z = 0.5,0.45,0.4,0.35,0.2,0.1,0.0; the point of convergence of 
tracks refers to z = 0.5. Data points and the dotted line is from 
Vikhlinin et al. (2002). 

cooling is dominant in the central region, as in cooling flow 
clusters, the direction of heating due to thermal conduction 
may reverse from being outward to inward. The central re- 
gion would then heat instead of getting cooled as discussed 
earlier in the paper. As a matter of fact, it has been specu- 
lated that radiative cooling may change the structure of the 
magnetic field and decrease the suppression of thermal con- 
duction in the central region (Soker & Sarazin 1990). It is 
also possible that both AGN heating and radiative cooling 
are important in cooling flow clusters; Ruszkowski & Begel- 
man (2002) found that heating due to thermal conduction 
suppressed only to the extent oi rj ^ 0.1, combined with 
AGN heating, may explain the recent observations of cool- 
ing flow clusters. 

We have also implicitly assumed the suppression factor 
to be homogenous throughout the intracluster gas, which 
may not always be valid, especially when cooling is impor- 
tant (sec above). Fabian, Voigt and Morris (2002) have dis- 
cussed the possibility that conduction is stronger in the cen- 
tral region than in the outer radii owing to the change in the 
structure of magnetic fields brought about by cooling in the 
centre. 



6 SUMMARY 

We have determined the time evolution of the temperature 
profile as a result of thermal conduction in the limit that 
the density profile is in steady state, beginning with tem- 
perature profiles obtained from numerical simulations (with 
thermal conduction). We have constrained the suppression 
factor for the thermal conduction to be of order < 10"^ 
by comparing the final temparature profile (after 10^" yr) 
with the observed profiles. We have also calculated the cor- 
responding evolution of the X-ray luminosity of clusters at 



a given temperature, but found the observed data too scat- 
tered to allow any useful constraint on thermal conduction. 
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